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Abstract 

We present a method for detecting neutrons using scintillating lithium gadolinium borate crystal 
shards in a plastic matrix while maintaining high gamma rejection. We have procured two cylin- 
drical detectors, 5"x5", containing 1% crystal by mass. Crystal shards have a typical dimension 
of 1 mm. One detector was made with scintillating plastic, and one with nonscintillating plastic. 
Pulse shape analysis was used to reject gamma ray backgrounds. The scintillating detector was 
measured to have an intrinsic fast fission neutron efficiency of 0.4% and a gamma sensitivity of less 
than 2.3 x 10~^, while the nonscintillating detector had a neutron efficiency of 0.7% and gamma 
sensitivity of (4.75 it 3.94) x 10^^. We determine that increasing the neutron detection efficiency by 
a factor of 2 will make the detector competitive with moderated ^He tubes, and we discuss several 
simple and straightforward methods for obtaining or surpassing such an improvement. We end with 
a discussion of possible applications, both for the scintillating-plastic and nonscintillating-plastic 
detectors. 
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1 I. INTRODUCTION 



2 Neutron detectors are useful in a variety of situations, from monitoring nuclear power 

3 plants |T], to imaging |2], to fissile material localization |3J, and so on. Depending on 

4 the specifics of the detector, they can be optimized for spectroscopy, directionality, scalar 

5 counting with or without timestamps, or some combination of the three. 

6 One detector technology that can be used for either counting or spectroscopy utilizes 

7 embedded inorganic scintillator crystals in a plastic matrix. This method was first proposed 

8 almost 10 years ago by Czirr et al. for the purposes of capture-gated neutron spectroscopy [1] . 

9 In that work, the crystal Czirr et al. used was lithium gadolinium borate (LGB) [5], with 

10 the chemical composition Li6Gd(B03)3, where the lithium and/or boron was enriched in the 

11 high-capture-cross-section isotopes ^Li and for various measurements. Capture-gated 

12 spectroscopy works by looking for a double light pulse, the first pulse coming from neutron 

13 recoils in the organic scintillator, and the second coming from the capture of that same 

14 neutron. The the case of LGB/plastic composite scintillators, interactions were differentiated 

15 by exploiting the very different scintillation decays times: a few ns for plastic, and 270 ns 

16 for the crystal. 

17 While it was possible for neutrons to recoil in the crystal rather than the plastic, this was 

18 unlikely in Ref. [1] because the detector was 90% plastic by mass. In a similar vein, it was 

19 also possible for the neutrons to capture on a proton in the plastic rather than the crystal, 

20 although this effect was subdominant because of the high macroscopic neutron capture cross- 

21 section of the crystal. A simple, though not absolutely precise, correlation could therefore be 

22 made, where short pulses were caused by neutron recoils, and long pulses caused by neutron 

23 captures. This ansatz extends to gamma recoils being typified by short pulses, similar to 

24 neutron recoils. 

25 Further particle discrimination is available using energy discrimination. For example, a 

26 neutron capture on ^Li liberates 4.78 MeV of kinetic energy. After scintillation quenching, 

27 this energy deposition in LGB produces the same number of optical photons as a 2.2 MeV 

28 electron This electron-equivalent energy is above most background radiation, with the 

29 ubiquitous 2.6 MeV gamma ray from ^os^l being a notable exception. Though the number of 

30 scintillation photons may be the same, a 2.2-MeV electron has a path length of approximately 

31 4.4 mm [6], so if the LGB crystals have a typical dimension of, for example, 1 mm, a 2.2-MeV 
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1 electron will create only ~500 keV of visible energy from the LGB shard. By comparison, the 

2 recoiling tritium resulting from a capture on ^Li has a range of only ~30 /im [7\ . Combining 

3 both energy analysis and pulse shape analysis is therefore capable of providing very strong 

4 discrimination between neutron captures and neutron/gamma recoils, as discussed below in 

5 Section HIl 

6 Additional work has been performed on LGB-embedded, scintillating-plastic detectors. 

7 Menaa et al. have evaluated 2" x 2" cylindrical LGB / scintillating-plastic detectors for 

8 possible use as a hand-held neutron spectometer [S]. Flaska, Pozzi, and Czirr have studied a 

9 5" diameter, 4" long detector from the standpoint of pulse shape discrimination [9] . Nelson 

10 and Bowden have studied such detectors to determine their suitability for anti-neutrino 

11 detection ITOl. 



12 The current work builds on the previous efforts in a number of ways. In Section II B 



13 we discuss methods to improve the pulse shape discrimination using additional information 

14 obtained from a digitized signal. In Section |III| we use Monte Carlo studies in an effort to 

15 increase neutron sensitivity while keeping gamma rejection at acceptable levels. Finally, in 



16 Section IV we briefly cover additional methods to improve detector performance. 



17 II. THE EXPERIMENTAL PROGRAM 

18 We began our studies with the same detector, data acquisition, and analysis as outlined in 

19 Ref. [in|. In that work, the detector, manufactured by MSI/Photogenics, was a cylindrical, 

20 5" x5" detector and 1% by mass LGB. The lithium and boron in the crystals were enriched to 

21 ~95% ®Li and ^°B, while the gadolinium isotopes were of natural abundances. The relevant 

22 nuclear reactions were therefore 

BR = 7% 

BR = 93% (1) 



plastic used was EJ-290, a 
5" Adit model B133D01W 



''Li + n -> 3h (2.75 MeV) + ^He (2.05 MeV) 
lOB + n ^ ^Li (1.0 MeV) + ^He (1.8 MeV) 

^ *^Li (0.83 MeV) + ^He (1.47 MeV) 
*^Li ^Li + 7 (0.48 MeV) 
^Gd + n ^ ^+lGd + 7 (~ 8 MeV cascade) 

23 The LGB shards had a typical dimension of 1 mm, and the 

24 scintillating polyvinyltoluene from Eljen Technology, Inc. Two 
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1 photomultiplier tubes, on either side of the detector, were connected to the data acquisition 

2 system. Attenuation in the hght signal was taken into account using standard analysis 

3 techniques, as described in Eq.(l) in Ref. [TU] . 

4 An analog data acquisition system was used to identify captures on ^Li by plotting the 

5 detector responses in a two-dimensional histogram, with the pulse shape parameter along 

6 the vertical axis and energy along the horizontal axis. Fig. [T] shows a clear signal island 

7 resulting from the capture of neutrons on ^Li. 

8 With this plot as a starting basis, we undertook several approaches to increasing the 

9 fidelity of the neutron capture signal, while increasing gamma rejection. We then mea- 

10 sured neutron efficiency and gamma rejection using both the scintillating-plastic and 

11 nonscintillating-plastic detectors. This work is detailed in the rest of this section. 
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FIG. 1. Results from the scintillating plastic detec- 
tor using an analog data acquisition system. This 
plot, along with an explanation of the pulse shape 
parameter, are from Ref. HQl- The captures in ^Li 
are clearly distinguishable as the island between 2 
and 3 MeVee and 0.8 to 1.0 on the vertical axis. 
A second island, much closer to the gamma back- 
ground, is also visible between 0.5 and 1 MeVee and 
0.6 to 0.8 on the vertical axis. This second island 
is associated with captures on ^^B. 



4 



1 A. Data acquisition 



2 In Ref. [in], an analog data acquisition system was used, allowing for a relatively simple 

3 pulse shape analysis using a ratio of the integral of the tail of the pulse to the integral of the 

4 full pulse. This tail / full ratio is the vertical axis of Fig. [T} Using this simple discrimination, 

5 there are several event topologies that can fake a neutron signal. For example, a low-energy 

6 electron recoil in the plastic may trigger the data acquisition system, and if an accidental 

7 high-energy electron recoil in the plastic falls within the tail gate time, the ratio of the tail 

8 integral to the / full integral may be close to 1, mimicking the tail / full ratio of a neutron 

9 capture on ^Li. 

10 To address this and other false neutron signals, we moved to a digital data acquisition sys- 

11 tem to provide finer control over the pulse shape analysis. We connected the photomultiplier 

12 tubes of both the scintillating-plastic and nonscintillating-plastic detectors to a Struck 3320 

13 digitizer, sampling at 200 MHz and 12 bits. The interface to the VME crate was a Struck 

14 3150 controller, and the data acquisition computer was running Linux. The data acquisition 

15 software was custom developed at Lawrence Livermore National Laboratory. Fig. [2] shows 

16 a photograph of the detector setup. 




FIG. 2. A photograph of the detector setup used in 
this work. Each detector had two photomultiplier 
tubes, which were run into an oscilloscope for test- 
ing purposes, and then on to the data acquisition 
system (not visible in this photo). The lead bricks 
between the detectors was used to reduce gamma 
cross-talk between the detectors. 
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1 The Struck 3320 firmware used allows for defining eight integral gates, four of which can 

2 be up to 512 samples each, and four of which can be up to 16 samples each. At a sample 

3 rate of 200 MHz, this translates to maximum gate times of 2560 ns for each gate in the first 

4 group of four, and 80 ns for each gate in the second group of four. The actual gate lengths 

5 used in the analysis were 50 samples for gates 1-4 and 16 samples for gates 5-8. Note that 

6 the start times of the gates can be arbitrarily set, allowing not only for gaps between gates, 

7 but overlapping gates if desired. Fig. [3] shows an idealized pulse from a neutron capture on 

8 ^Li, and how the gates were set up to subdivide the pulse. For every event, we used the 

9 baseline values as measured by gates 7 and 8 to obtain baseline-subtracted integrals on gates 

10 1-6. This baseline subtraction of course took into account the differing gate widths. 

11 The best possible pulse shape discrimination would have been possible if we recorded 

12 each entire pulse and ran pulse shape analysis on that, rather than relying on simple gate 

13 integrals. This would have introduced a great amount of computing overhead, though, as 

14 400 points would have to be transferred to the DAQ computer and analyzed, rather than 
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FIG. 3. Setup of the 8 integral gates on the Struck 
3320. Gates 7 and 8 measure the baseline, gates 5 
and 6 measure the rising edge, and gates 1-4 are 
used for pulse shape analysis. Summing gates 1-6 
and subtracting the baseline contribution gives the 
full integral of the pulse. The curve is an idealized 
signal from the LGB crystal for clarity of explana- 
tion, and is not an experimental pulse. 
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1 just 8. Because we used the gate integrals, we were able to maintain DAQ livetime above 

2 99%, even when the event rate was on the order of 5 kHz. 



3 B. Improving the pulse shape analysis 

4 Once the digital data acquisition was set up, we re-created the analog pulse shape analysis 

5 by constructing a tail / full ratio: 

Integral of gates 1-4 

Analog-like pulse shape parameter = — [2) 

Integral of gates 1-6 

6 The results of this analysis for both detectors is shown in Fig. |4j The data in these plots 

7 were obtained using a bare, 75.4 kBq 252Qf source (8690 neutrons/second) suspended ap- 

8 proximately 3 feet above the detectors. Data was accumulated over 20 hours, and the live 

9 time was over 99%. For both detectors the ^Li island is clearly visible. The very broad 

10 gamma band between pulse shape parameter 0.1 and 0.2 in the scintillating-plastic detector 

11 is absent from the nonscintillating-plastic detector, while the amplitudes of the ^Li islands 

12 are of similar size. 

13 We have identified multiple event topologies that lead to differing effects when construct- 

14 ing the two-dimensional plots such as are shown in Fig. |4j Fig. |5] shows possible pulses 

15 from the scintillating-plastic LGB detector. The ^°B capture shape results from the cascade 

16 gamma depositing energy in the plastic while the ion recoils remain within the LGB crystal. 

17 This pulse topology can also result if a neutron capture within an LGB shard results in 

18 ions escaping the crystal, such that energy is deposited in both the plastic and the crystal. 

19 Similarly, a gamma recoil can result in an electron depositing energy in both the plastic and 

20 a crystal shard. 
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Pulse shape parameter vs. Energy for the scintillating detector 




Energy (arb. units) 



Pulse shape parameter vs. Energy for the Nonscintillating detector 




Energy (arb. units) 



FIG. 4. Analog-like data analysis of the 
scintillating-plastic (top) and nonscintillating- 
plastic (bottom) detectors. The tail gate was not 
identical to that used in the previous work, so the 
pulse shape parameter calculations will not be iden- 
tical. Compare the top plot to Fig. [TJ The analog- 
like data analysis from this work differs from that 
shown in Fig. [T] because of differences in PMT re- 
sponse, integral gates, and threshold levels. 
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FIG. 5. Possible pulses from the scintillating- 
plastic LGB detector. These pulses are idealized 
pulses, and not from experimental data. See text 
for details. 
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1 Fig. |5] also demonstrates the maximum data rate of these composite detectors. The 

2 Struck 3320 has data buffers that allow it to digitize a signal while transferring previously- 

3 accumulated data to a computer. This results in very low dead time, as discussed in Sec- 



4 tion lie Because each pulse is digitized for ~1 jjs, however, it is possible for multiple events 

5 to arrive within a single digitization time window. For example, the second gamma pulse in 

6 the bottom curve of Fig. [5] could just as well be a neutron-capture event. If this were the 

7 case, the second pulse would be lost as part of event pile up. Because of the /xs-scale digi- 

8 tization window, if the event rate were greater than ~100 kHz, pile up would significantly 

9 increase. 

10 Using the simple pulse shape analysis, an accidental gamma coincidence can create a 

11 background to the neutron capture signal. Referring to the bottom curve of Fig. [5| we note 

12 that the "tail" of the pulse, integrated from 500 ns to the end of the pulse, is a sizeable 

13 proportion of the full integral. This has the tendency to pull the pulse shape parameter 

14 closer to unity, and the event ends up encroaching on the neutron capture regions. We can 

15 avoid these kinds of backgrounds by making additional use of the gate integrals depicted in 

le Fig.|3} 

17 The approach we used to enhance the pulse shape discrimination was to use the ratios of 



(Gate 5 + Gate 6) / Gate 1 
Gate 1 / Gate 2 

(3) 

Gate 2 / Gate 3 
Gate 3 / Gate 4 

18 We measured these ratios using events from the clearly-defined ®Li-capture islands from 

19 Fig. |4] for both the scintillating-plastic and nonscintillating-plastic detectors. Fig. [6] shows 

20 the distribution of gate ratios as measured using the scintillating-plastic data. We then 

21 re-analyzed the full data sets, requiring all events have gate ratios as per the ^Li-capture 

22 islands. We then plotted the remaining events on a two-dimensional histogram analogous 

23 to those in Fig. |4]-with the simple pulse shape parameter on the vertical axis-for direct 

24 comparison between the simple and updated pulse shape discrimination procedures. These 

25 plots are shown in Fig. [7} The success of our approach emphasizes the requirement that the 

26 LGB have a sufficiently long scintillation decay time. 

27 Even a cursory comparison of Fig. |4] and Fig. [7| shows the efficacy with which the updated 
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1 pulse shape discrimination removes background events. 



Ratios of scintillating gate integrals 




FIG. 6. Gate ratios for the scintihating-plastic de- 
tector. The curves' asymmetry prevented us from 
fitting Gaussian curves to the data. We incorpo- 
rated manual ratio limits, however, that included 
at least 90% of the curves' areas. For example, for 
the (Gate 5 -|- Gate 6) / Gate 1 ratio, we used 
limits of 1.5 and 2.2 (shown by the vertical bars), 
incorporating 92% of the available area. Analogous 
limits were identified for the other three ratios. 
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Pulse shape parameter vs. Energy for the scintillating detector 




Energy (arb. units) 



Pulse shape parameter vs. Energy for the Nonscintillating detector 




Energy (arb. units) 



FIG. 7. Data after applying the updated pulse 
shape discrimination procedure. Compare to 
Fig. |4j The ^Li capture regions are unattenuated, 
while the gamma backgrounds have been drasti- 
cally reduced. Note the second island in the lower 
plot, which is associated with captures on ^"^B. This 
island does not appear in the upper plot because of 
the 0.48 MeV gamma that accompanies ^*^B cap- 
tures. 
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C. Measurements of neutron and gamma sensitivity 



We measured the intrinsic neutron capture efficiency and gamma sensitivity of our two 



detectors by acquiring three data sets. The ffist was the data set described in Section |II B 
a 20-hour data set with a bare 252Qf source approximately 3 feet from the detectors. The 
second data set was a 20-hour data set with a 14.0-kBq ^^^Th source, also suspended ap- 
proximately 3 feet from the detectors. The final data set was a 20-hour background data 
set with no sources present. All three data sets had live times over 99%. The maximum 
trigger rate was less than 10 kHz, ensuring that event pileup was not a concern. Given the 
source activities, length of data accumulation, and source distance from the detectors, we 
calculated a total neutron flux of 9.60 x 10^ from the 252q£ source and a 2614-keV gamma 
flux of 1.39 X 10^ from the ^^^Th source. 

We would like to comment on the use of a ^^^Th to measure the effect of gamma rays on our 
composite detectors. Typically, a ^'^^Cs source would be used to make such a measurement. 
Given the 2.2-MeV energy equivalent from a neutron capture on ^Li, the 662-keV gamma 
rays from ^^^Cs would virtually never create a background to the neutron capture signal, 
even if pulse shape discrimination were not used. Three gamma rays would have to interact 
within the LGB shards simultaneously, and deposit nearly 100% of their energy within the 
shards. The probability of this is vanishingly small, even in a high-rate environment. To have 
a reasonable expectation of any gamma background at all, we needed to use a source with 
gamma rays at least 2.2 MeV. Of course, the gamma rays resulting from neutron captures 
on hydrogen are one possible source, but we wanted to separate the neutron and gamma 
signals. Thus we made use of the ^^^Th source, which contains within its decay chain ^°^T1, 
and its 2614-keV gamma ray. 

In the end, our gamma ray sensitivity measurements were more stringent than typical 
gamma sensitivity measurements performed on, for example, ^He proportional counter tubes. 

For all three data sets we applied the updated pulse shape discrimination, then utilized an 
additional cut in the simple tail / full ratio of between 0.4 and 0.6 for the scintillating-plastic 
detector, and between 0.4 and 0.7 for the nonscintillating-plastic detector. We projected the 
resulting data onto the x-axis and performed a bin-by-bin subtraction of the background 
histogram from both the gamma and neutron histograms. The results are shown in Fig. [8] 

Normalizing the integrals in the neutron signal regions by the total flux of either neutrons 
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Scintillating-plastic detector spectra (Background subtracted) 
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FIG. 8. Background-subtracted data for neutron 
and gamma sources. We fit Gaussian curves to the 
neutron data, and integrated over ±3a for both 
the neutron and gamma curves. The integrals are 
shown on the plots. Because there was no normal- 
ization applied to the curves before subtracting, the 
integrals summed to integer counts. For both plots, 
the dashed vertical lines show the region of integra- 
tion. 



1 or 2614-keV gammas, we calculated the intrinsic neutron efficiency and gamma sensitivity 

2 of these detectors. Results are shown in Table [B 

3 We can compare the intrinsic efficiency of our LGB composite detectors to fast neutrons 
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TABLE I. Intrinsic gamma and fission neutron sensitivity of tlie two detectors. Note that because 
the gamma curve in the scintillating-plastic detector integrates to less than 0, we can only apply a 
la upper limit on gamma sensitivity. 



Detector 


Neutron sensitivity 


Gamma sensitivity 


Scintillating-plastic 
Nonscintillating-plastic 


(0.416 ± 0.007)% 
(0.683 ± 0.009)% 


< 2.3 X 10-9 
(4.75 ± 3.94) X 10-^ 



with that of moderated ^He tubes. East and Walton constructed an array of five ^He 
tubes [H]. Each tube was 8" long, 1.5" in diameter, and pressurized to 6 atmospheres. 
The tubes were embedded in a polyethylene moderator 9" in diameter, as shown in Fig. |9| 
The five tubes were fed to a single preamplifier, and the total intrinsic efficiency for fast 
neutrons was measured to be (11.5 ±0.5)%. Our 5" detectors would fit very cleanly into the 
9" cylinder, providing for 2" of polyethylene moderator. Menaa et al. measured a 10-fold 
increase in fast neutron efficiency for their LGB composite detector simply by using 2" of 
polyethylene moderator [8]. We would therefore expect our nonscintillating-plastic detector 
to have an intrinsic efficiency for fast fission neutrons of ~7%. Increasing the performance 
of our detectors by a factor of two would make them competitive with ^He-based systems. 
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INCH OIA X 3.5 INCH DEEP 
HOLES fTWELVE) DRILLED 
INTO POLYETHYLENE CYLINDER 



FIG. 9. Moderated ^He tube array for fast neutron 
detection. East and Walton measured an intrinsic 
fast neutron efficiency to be (11.5ib0.5)%. Tlie cen- 
tral active volume is approximately 5" in diameter, 
the same as our composite LGB detectors, allowing 
for easy comparison between the detectors. Figure 
taken from Ref. jllj. 
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III. OPTIMIZING CRYSTAL SIZE AND CONTENT 



The detectors utilized in this work were intended to demonstrate the feasibihty of com- 
posite detectors, and as such they were not optimized for neutron sensitivity or gamma 
rejection. In this section we discuss possible optimizations to increase neutron sensitivity 
while maintaining or reducing gamma sensitivity. 



A. Crystal size 



We begin the optimization analysis from the point of view of crystal self-shielding. The 
LGB crystals provide a signal when a neutron captures on the ^Li nuclei, and the capture 



depth of the neutrons is shown in Fig. [TOj At thermal energies, the capture depth is on 
the order of 10-20 fim, although neutrons are not necessarily completely thermalized by the 
time they capture. Regardless, the capture depth will tend to be small compared to the 
typical shard dimension of 1 mm, so the shards will exhibit considerable self-shielding. The 
captures occur primarily in the outer edge of the crystals, and the inner volumes are more 
or less unused. 

From the point of view of self-shielding, the smaller the crystal shards, the more effi- 
ciently the LGB crystals capture neutrons, even with the same total mass. Recall, however, 
that the recoiling ions have a range of approximately 30 fim, so if the crystals become too 




10'^ 10" 10'' 10"=^ 10"^ 1 10 10^ 10^ 10" 10^ 10^ 10^ 

Neutron energy (eV) 



FIG. 10. Attenuation depth for neutrons in the 
LGB crystal. This crystal, as described above, is 
enriched to 95% in both ^Li and ^'^B. 



17 



small, we again potentially lose efficiency because less scintillation light will be produced in 
the shard. There are therefore multiple parameters that need to be considered when calcu- 
lating the optimal shard dimension, including the primary neutron spectrum, efficiency of 
thermalization, capture cross-section of both the plastic and the LGB crystal, angle of ion 
recoil with respect to the shard surface, and so on. These multiple related variables make 
an analytic calculation intractable, so we turn to a Monte Carlo calculation to determine 
the optimal shard dimension. 

To determine an appropriate range over which to run the Monte Carlo calculations should 
be run, we perform a simple one-dimensional calculation to predict the optimal dimension. 
We assumed a constant-mass LGB crystal with a uniform flux of monoenergetic neutrons 



normal to one face of the crystal (Fig. 11 ). D is the depth of the crystals along the direction 
of the incident neutron flux. 

We define A as the capture depth of a neutron, which in this treatment is taken to be 
a constant, in accordance with the simplification of monoenergetic neutrons. The capture 
probability for a single neutron entering the crystal is therefore 1 — e~^f^. The total number 
of neutrons entering the crystal is the neutron flux multiplied by the cross-sectional area A. 
The smaller the value for D, the larger the total flux. The total capture efficiency therefore 



Neutrons 




Capture location x 



FIG. 11. Setup for calculating a simple, one- 
dimensional optimization of the LGB shard size. 
D is the length of the crystal, and A is its cross- 
sectional area. Thus D and A are inversely propor- 
tional. R is the range of the recoil ions. Because 
this is a one-dimensional calculation, the angle of 
ion recoil is irrelevant. 
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goes as 



Efficiency ~ (l - e^/^) /D (4) 

Next, we note that if the capture location within the crystal is greater than a distance 
R from either end of the crystal, the recoiling ions will deposit 100% of their energy in the 
crystal. Similarly, ii D < 2R, the ions will never deposit full energy in the crystals. We 
therefore integrate Eq. Q over the capture location, x, from to D to obtain the probability 
of the recoiling ions depositing their full energy in the crystals, assuming a minimum value 
of 2R for D. Assuming conservative values of 40 fim for R and 200 fim for A, the efficiency 



curve is shown in Fig. 12 Within the limits of this simple approximation, we expect the 



optimal shard size to be on the order of hundreds of microns. 



B. Total crystal content 

A second optimization to be performed using Monte Carlo calculations is the total crystal 
content of the composite detector. The detectors used in this study were 1% by mass LGB, 
and increasing crystal content would increase neutron sensitivity. Increasing crystal content, 
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FIG. 12. Simple, one-dimensional optimization of 
the crystal shard size for i? = 40 ^m and A = 
200 ^m. The maximum occurs near 300 fj,m, so our 
Monte Carlo calculations will approximately span 
the range 100 fim to 1 mm. 
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1 however, might also increase gamma sensitivity, and thus the Monte Carlo calculations must 

2 be run with both fission neutrons and gamma rays. 

3 The effects of total crystal content on both neutron and gamma sensitivity can be some- 

4 what subtle. If the crystal content is low, the neutron sensitivity should be expected to 

5 scale linearly with it. The gamma sensitivity, however, is a strong function not only of total 

6 crystal content, but crystal dimension as well. Recall from Section |T] that the range of a 

7 2.2-MeV electron in the LGB crystal is 4.4 mm, while the range of the tritium ion from 

8 capture on ^Li is only about 30 /im. This implies that if the shards are small, the gamma 

9 background to the neutron capture signal can be kept small. Unfortunately, if the crystal 

10 content is sufficiently large, it is possible for a recoiling electron to traverse multiple LGB 

11 shards. In this case, even if each shard emits only a fraction of the light equivalent of a 2.2 

12 MeV of electron energy, the sum total of the light emitted by all the shards may create a 

13 background to the neutron signal. 

14 C. Monte Carlo Results 

15 For our Monte Carlo calculations, we used LUXSim, a simulations package based on 

16 Geant4 |T2]. Geant4 is a C++-based framework for simulating the effects of various kinds of 

17 radiation in different materials, driven by both physics models as well as data [I3]. Geant4 

18 is capable of handling electrons, gammas, neutrons, and ions. 

19 For the calculations, we used Geant4, version 4.9.3.p01. Unfortunately, there was also a 

20 documented bug in Geant4 that results in incorrect ion recoil energies resulting from neutron 

21 captures on ®Li [13] . For our simulations, we implemented a correction so that the recoiling 

22 tritium and helium nuclei had the correct energies. As a side note, there was also a bug for 

23 captures on ^°B, although since our analysis only relied on captures on ^Li, the effects of the 

24 ^°B bug could be safely ignored. 

25 The LGB crystal shards in our simulations were all placed individually inside the plastic 

26 cylinder. This was to ensure there was no overlap between the crystals, thereby guaranteeing 

27 the total crystal content is well-defined. Unfortunately, this resulted in an exceedingly large 

28 number of shards to place, each one having to be checked against overlap with any others. 

29 To make the placement more tractable, the cylinder was sub-divided into layers, and a 

30 proportionately smaller number of crystal shards were placed in the layer. The layer itself 
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was then duplicated to create the full cylinder. An image of such a construction is shown in 
FiK.ITSl 



For each simulation data set, we recorded the total energy deposition in the plastic and the 
crystal, as well as the particle that created that energy deposition. If the particles depositing 
the energy came from neutron captures on ^Li, their energy deposition was quenched by a 
factor of 0.46 (the ratio of 2.2 MeV / 4.78 MeV discussed in Section [l]). A finite resolution of 
9% was applied to the energy depositions, to match the width of the Gaussian peak in the 
upper plot of Fig. [sj Finally, a 3a range around the centroid of the ^Li peak was defined, 
and any counts within this range was tabulated. An example of this entire process is shown 



in Fig. 14, which shows the results from a simulation with 1-mm crystal shards, and 1% 
LGB content, with a fission neutron source. 

We ran many simulations to fill out the parameter space of interest. We varied the 
primary particles between 252Qf fission neutrons, and ^os^l decays to provide the 2614-keV 
gamma rays from the ^^^Th decay chain. We varied the shard dimension from 100 fim to 
1 mm in 100-/im steps, and we varied the total mass of the LGB crystals from 1% to 10% 



by mass, in 1% steps. The results for the fission neutron source are shown in Fig. 15 In 




FIG. 13. Example Geantl geometry of the com- 
posite detector. In this example, the shard size is 
0.5 mm, the total content is 1%, and there are 5 
layers. Note that the shards, while having identi- 
cal dimensions, are placed in the cylinder in ran- 
dom orientation. The layers themselves, though, 
have identical placement of the shards. 
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Energy (MeV) 

FIG. 14. Sample plot of analysis of a Geant4 sim- 
ulation with a neutron source. The unquenched 
energy depositions range up to 10-12 MeV because 
of the sum of the energy of the primary neutron and 
the energy of the 8-MeV gamma cascade resulting 
from neutrons capturing on gadolinium nuclei. See 
text for details on each curve. 



this figure, we see that while the curves are much flatter than the simple curve from Fig. [12 
the optimal shard size lies generally in the 200-/im to 600-/im range. We also note that, 
based on these Monte Carlo simulations, we can increase the neutron sensitivity by roughly a 
factor of 5 by using 10% LGB content and 500-/xm-sized shards. We note also, though, that 
given the somewhat fiat response, the intrinsic neutron efficiency is somewhat insensitive to 
precise shard dimensions, provided the shards themselves are larger than ~300 fim. 

The question remains, however, of the effect the crystal size and content on the gamma 
sensitivity. Using the ^os^jj (;iata, we counted the number of events in the ^Li neutron capture 



region. Sample analysis plots similar to that in Fig. 14 are shown in Fig. 16 They show 
that, even with the same total LGB content, the shard size can still have a substantial 
effect on the potential gamma sensitivity to these composite detectors. The final plot from 



our gamma sensitivity simulations is shown in Fig. 17 From this plot we conclude that at 
the optimal shard size for neutron captures, 500 /im, the gamma sensitivity is consistently 
reduced relative to a shard size of 1 mm. This gives us confidence that the LGB content 
can be increased to the 10% range without increasing gamma sensitivity. 
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FIG. 15. Optimizing the shard size and total LGB 
content for the composite detector. Each curve rep- 
resents the amount of LGB in the detector by mass. 
The detectors used in this study had a shard size of 
~1 mm, and a content of 1%. Computing hmita- 
tions prevented us from running simulations with 
100-;um shard size above 4% content. 

One final note on crystal optimization is the effects of small shards on optical trans- 
parency. The LGB has an index of refraction of 1.66, compared to the EJ-290's value of 
1.58. Because of this refractive index mismatch, a plastic matrix with a high LGB con- 
tent may be optically cloudy (see, e.g.. Fig. 1 from Ref. [8], where the detectors were 27% 
LGB by mass). A full treatment of the optimization of LGB content would have to include 
scintillation photon attenuation. This attenuation must be measured for any given final 
configuration, and extrapolations or interpolations made from there to determine optimal 
LGB content given the overall dimensions of the plastic matrix. 
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FIG. 16. Sample analysis plots of the Geant4 sim- 
ulations with the 2*^^X1 source. For both of these 
plots, the total LGB content was 10%. For the top 
plot, however, the shard size was 0.2 mm, and for 
the bottom plot the shard size was 1 mm. We see 
that the shard size can greatly affect the number of 
gamma counts in the ®Li neutron capture region. 
Indeed, the top plot had zero counts in this region, 
while the bottom plot had 8 counts. The maximum 
possible total energy in the gamma rays from ^oS'pj 
decays is 4.4 MeV, although the branching ratios 
to those levels are very small. The branching ratio 
to the 3.96-MeV energy level is 3%, which is the 
approximate endpoint of the energy depositions in 
the plastic. 
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FIG. 17. Gamma sensitivity in the ^Li neutron 



capture region. This plot is analogous to Fig. 15 
but with a ^oS'pj source instead of a fission neutron 
source. There were no counts at all for LGB content 
lower than 4%. While the statistics are low, it is 
clear to see that, regardless of the LGB content, 
reducing the shard size from 1 mm to 0.5 mm will 
not increase the gamma sensitivity. 
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1 IV. ADDITIONAL IMPROVEMENTS 



2 In addition to optimizing the crystal content and shard size, we have identified other 

3 avenues of increasing detector performance. The goal is to increase neutron sensitivity while 

4 maintaining the very high level of gamma rejection. 

5 A. Changes to the data acquisition system 

6 One method is to change the gate sizes of the data acquisition software. This could be 

7 done through a firmware update, although the code itself would have to be contracted from 

8 Struck. Using the eight available gates, it would be beneficial to change their maximum 

9 sample sizes. One gate would be used to measure the baseline value on an event-by-event 

10 basis, and 50 samples would suffice. The LGB scintillation pulses only need to be digitized 

11 for 5 times the decay constant, or ~1350 ns. Subdividing this pulse by 7 gives roughly 39 

12 samples per gate. As it is, having 4 gates limited to 16 samples while 4 gates have 512 

13 samples available means there is a fair amount of flexibility in shuffling around the available 

14 space. Certainly, 8 gates with a maximum sample size of 50 would suffice for this application. 

15 With the 7 data points, one gate would cover the rising edge of the pulse. The other six 

16 samples could then be used to identify a very clean exponential decay. Indeed, to perform 

17 better pulse shape analysis, an exponential curve could be fit to the six decay data points, 

18 with the time constant a fixed parameter. The only free parameter in the fit would therefore 

19 be the amplitude of the pulse. The pulse shape discrimination would them come about via 

20 the reduced ^ of the fit, where an associated p-value between, for example, 0.05 and 0.95 

21 could be counted as a good fit. Anything else would be rejected as pileup, a gamma event, 

22 a muon event, or some other spurious background. 

23 An additional benefit of this fitting approach would be the reconstruction of the full pulse 

24 integral. Once the fit amplitude is known, the energy resolution could conceivably increase, 

25 leading to a smaller neutron capture region, and therefore reduced background. 

26 B. Alternatives to lithium gadolinium borate 

27 The analysis in this work relied entirely on neutron captures on ^Li within the LGB crys- 

28 tal. Recall that the lithium was enriched in ^Li and the boron in ^°B, while the gadolinium 
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1 was of natural isotopic content. In Ref. [10], the authors showed that using gadohnium and 

2 boron depleted in their high-capture-cross-section isotopes, the number of captures on ^Li 

3 could increase by a factor of 6. While depleted boron is not difficult to obtain, depleted 

4 gadolinium, unfortunately, is. 

5 We therefore identified materials other than LGB that could potentially be of use. These 

6 materials should have the following characteristics to make them suitable alternatives to 

7 LGB: 

8 • High lithium content 

9 • Long scintillation decay constant 

10 • No isotopes to compete with the lithium for neutron captures 

11 • High light output 

12 We have identified a few possible materials for future testing, shown in Table [11} Of these 

13 materials, the most directly compatible with the results from this work is LYB, since it 

14 merely replaces the gadolinium with yttrium. Given this simple extrapolation, we estimate 

15 the lithium capture rate would increase by a factor of 6 over LGB, in accordance with results 

16 from [To]. While the LYB light yield is small, this is not anticipated to be a strong issue, 

17 as the light yield on the scintillating plastic could be reduced a similar amount, while the 

18 gain on the PMTs could be turned up to compensate for the loss of light. If the refractive 

19 index of the plastic were matched to that of the inorganic crystal, the light collection would 

20 be further improved. 

21 Lithium glass presents another attractive alternative, for two reasons. One is that its 

22 index of refraction (n = 1.56), is even closer to that of the plastic (n = 1.58) than LGB 

23 {n = 1.66). This would allow for greater light propagation, counteracting the lower light 

24 output. Another reason is that none of its component materials would require enrichment, 

25 as neither oxygen nor silicon have isotopes with high neutron-capture cross-sections. 

26 LiF has the highest concentration of lithium atoms, but unfortunately its decay constant 

27 is far too long to be a viable LGB replacement. It may be possible that a dopant other 

28 than tungsten would result in shorter scintillation decay time, but such a combination is 

29 unknown to the authors. 

30 CLYB is a known scintillator that has intrinsic neutron / gamma discrimination proper- 

31 ties. Unfortunately, its lithium content is low (1 out of 10 atoms) relative to LGB (6 out of 
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TABLE II. Possible alternatives to LGB. 



Material name 


Composition 


Light Output 
rel. to LGB 


Scint. 

r (ns) 


Reference 


LYB 


Li6Y(B03)3 


0.09 


100 


[15J 


Li glass 


Variable 


0.29 


62 


il6j 


LiF 


LiF(W) 


0.14 


40000 


[17J 


CLYB 


CsaLiYBre 


2.21 


85 


[18J 


LTO 


Li3Ta04 


2.14 


6100 


[19J 


LTB 


LiBaOs 


0.05 


< 10 


[20J 



19 atoms). Still, its high light output and long scintillation decay constant may make it a 
worthwhile LGB replacement . LTO has a high lithium content and high light output, but 
its scintillation decay time is too long to make it of use in a high-rate environment, as any 
detector using LTO would be limited to neutron capture rate of roughly 30 kHz. The very 
short decay constant of LTB would make its signal near indistinguishable from that of the 
plastic scintillator, or even Cherenkov light produced in the nonscintillating plastic detector. 



V. APPLICATIONS 



Several applications are available for the detectors discussed in this work. One of the 
more immediate applications would be as a neutron scalar detector. Given that the contin- 
uing shortage in the ^He supply is not expected to diminish in the future, these detectors. 



with the improvements discussed in Sections III and IV B could be used as a replacement 
for ^He-based proportional counters. These detectors would have intrinsic efficiency for de- 
tecting fission neutrons in the 10-20% range, while maintaining a strong gamma rejection 
capability of at least 10^^. Either the scintillating-plastic or nonscintillating-plastic detec- 
tor would work for this application, although for a high gamma radiation environment, the 
nonscintillating-plastic detector would be preferable to reduce the overall trigger rate. 

Another application would be an improved capture-gated spectrometer, as described 
in [4j. The higher neutron capture rate afforded the various optimizations, combined with 



the stronger pulse shape analysis described in Section [IVA[ would allow for a more accurate, 
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1 and higher-efficiency, spectrometer. Such a detector would be made of the scintillating 

2 plastic. 

3 Finally, the scintillating-plastic detector could possibly be used for anti-neutrino detec- 

4 tion. These detectors rely on a hydrogenous material to provide a target for inverse beta 

5 decays, and then look for a double signal coming from the prompt positron annihilation, 

6 followed by the neutron capture [21]. Previous anti-neutrino detectors allow for some event 

7 selection based on energy deposition, but they largely do not utilize pulse shape in differen- 

8 tiating between neutron captures and gamma recoils. Using a scintillating-plastic detector 

9 of the kind described in this work would reduce background from both accidental gamma 

10 coincidences as well as multiple-neutron capture resulting from a single cosmic ray spallation 

11 event near the detector. 



12 VI. CONCLUSIONS 

13 Lithium gadolinium borate crystals do not have any intrinsic neutron / gamma discrim- 

14 ination, but if small shards are embedded in a plastic matrix, the resulting signal can be 

15 used to differentiate neutron captures and gamma recoils with a high degree of confidence. 

16 The discrimination relies on accurate pulse shape identification, specifically between the 

17 fast-decaying signal associated with gamma recoils and the slow-decaying signal associated 

18 with neutron capture. Energy analysis also provides for particle discrimination, leading to 

19 gamma sensitivity on the order of 10~^. 

20 The intrinsic neutron efficiency of the detectors studied in this work is on the order 

21 of 0.4-0.7%, although they were not optimized for neutron detection. Four avenues exist 



22 for increasing neutron sensitivity, outlined in Table III These efficiencies are of course not 

23 orthogonal, although we expect that if all methods were employed, the performance of these 

24 detectors would surpass that of ^He-based detectors. From a quick sample survey, the most 

25 promising LGB alternatives are lithium yttrium borate and lithium glass. The gamma 

26 sensitivity could be kept low by decreasing the shard size, as well as using a more finely- 

27 tuned data acquisition system designed for this kind of composite detector. Even a factor 

28 of 15 improvement on the current neutron detection sensitivity would give these detectors a 

29 neutron efficiency comparable to moderated ^He tubes. 

30 Apphcations for these detectors include neutron counting (specifically ^He-based detector 
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TABLE III. Optimizations to the LGB detectors, and their anticipated increase in fast neutron 
efficiency over the baseUne performance. 



Method 


Factor improvement over baseline 


2" polyethylene moderator 


10 


Alternative crystal 


6 


Increased total crystal content 


4 


Smaller shards to limit self-shielding 


1.5 



replacement), capture-gated neutron spectroscopy, and anti-neutrino detection. 

This work performed under the auspices of the U.S. Department of Energy by Lawrence 
Livermore National Laboratory under Contract DE-AC52-07NA27344. LLNL-JRNL-499156. 
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